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ABSTRACT 

We have performed a detailed analysis of the local super-massive black-hole (SMBH) 
mass function based on both kinematic and photometric data and derived an accurate 
analytical fit in the range 10 6 M B n/M e < 5 x 10 9 . We find a total SMBH mass 
density of (4.2 ± 1.1) x 10 5 M Q /Mpc 3 , about 25% of which is contributed by SMBHs 
residing in bulges of late type galaxies. Exploiting up-to-date luminosity functions 
of hard X-ray and optically selected AGNs, we have studied the accretion history of 
the SMBH population. If most of the accretion happens at constant Mbh/-Mbh, as 
in the case of Eddington limited accretion and consistent with recent observational 
estimates, the local SMBH mass function is fully accounted for by mass accreted by X- 
ray selected AGNs, with bolometric corrections indicated by current observations and 
a standard mass-to-light conversion efficiency e ~ 10%. The analysis of the accretion 
history highlights that the most massive BHs (associated to bright optical QSOs) 
accreted their mass faster and at higher redshifts (typically at z > 1.5), while the 
lower mass BHs responsible for most of the hard X-ray background have mostly grown 
at z < 1.5. The accreted mass function matches the local SMBH mass function if, 
during the main accretion phases, e ~ 0.09 (+0.04, —0.03) and the Eddington ratio 
A = L/L^dd — 0.3 (+0.3, —0.1) (68% confidence errors). The visibility time, during 
which AGNs are luminous enough to be detected by the currently available X-ray 
surveys, ranges from ~ 0.1 Gyr for present day BH masses Mg H ~ 10 6 Mq to ~ 0.3 Gyr 
for Mg H ^ 10 9 Mq. The mass accreted during luminous phases is ^ 25-30% even if we 
assume extreme values of e (e ~ 0.3 — 0.4). An unlikely fine tuning of the parameters 
would be required to account for the local SMBH mass function accomodating a 
dominant contribution from 'dark' BH growth (due, e.g., to BH coalescence). 

Key words: black hole physics - galaxies: active - galaxies: nuclei - galaxies: evolu- 
tion - quasars: accretion - cosmology: miscellaneous 



1 INTRODUCTION 

The paradigm that quasars and, more generally, Active 
Galactic Nuclei (AGNs) are powered by mass accretion 
onto a super-massive black hole (SMBH) proposed long 
ago (Salpeter 1964; Zeldovich & Novikov 1969; Lynden-Bell 
1969) has got very strong support from spectroscopic and 
photometric studies of the stellar and gas dynamics in the 
very central regions of local spheroidal galaxies and promi- 
nent bulges. These studies established that in most, if not 
all, galaxies observed with high enough sensitivity a central 
massive dark object (MDO) is present with a well defined 
relationship between the MDO mass and the mass or the 
velocity dispersion of the host galaxy spheroidal component 
(Kormendy & Richstone 1995; Magorrian et al. 1998; Geb- 



hardt et al. 2000; Ferrarese & Merritt 2000; Tremaine et 
al. 2002; Kormendy 2003). Although there is no direct evi- 
dence that all MDOs are black holes (BHs) , the evidence for 
a singularity is actually very tight in the Galaxy (Schodel et 
al. 2002; Ghez et al. 2003) and alternative explanations are 
severely constrained in NGC 4258 (Miyoshi et al. 1995; see 
also e.g. Kormendy 2003). 

Soltan (1982) showed that, in the framework of the 
above paradigm, the total accreted mass density can be in- 
ferred from the observed QSO/AGN counts. The basic ingre- 
dients of the calculation are i) the bolometric correction fcbol, 
ii) the effective redshift and the corresponding K-correction, 
and iii) the mass to radiation conversion efficiency e. When 
more precise luminosity functions of QSO/AGN became 
available, Chokshi & Turner (1992) presented a first esti- 
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mate of the accreted mass density and derived constraints 
on the corresponding mass function (MF). Small & Bland- 
ford (1992) tried to relate the luminosity function (LF) of 
optically selected AGNs to the local SMBH MF in galaxies 
for several possible accretion histories. Salucci et al. (1999) 
showed that the LF of the galaxy spheroids (encompassing 
E and SO galaxies, and bulges of spiral galaxies) combined 
with the relationship between the spheroid and the central 
MDO mass allows an accurate evaluation of the local MF 
of the SMBHs. They concluded that the distribution of the 
mass that fuelled the nuclear activity, as traced by the LFs 
of QSOs and of AGNs contributing the X-ray background, 
matches the local SMBH MF, provided that e ~ 0.1 and the 
ratio of bolometric to Eddington luminosity A = Lboi/^Edd 
declines with luminosity and/or look-back time. It was also 
shown that the high mass tail of the local MF is due to the 
remnants of the bright optically selected QSOs, while at low 
masses the remnants of the AGNs producing most of the 
X-ray background largely dominate. 

Recent estimates of the local MF also exploited the rela- 
tionship between SMBH mass and velocity dispersion of the 
host spheroid (see e.g. Yu & Tremaine 2002; Aller & Rich- 
stone 2002; McLure & Dunlop 2003; Marconi et al. 2004), 
to reach, however, somewhat discrepant conclusions. 

The aim of this paper is to compare the local MF of 
SMBHs with the MF of the mass accreted during the nuclear 
activity (AMF), in order to shed light on important open 
questions such as: 

• is there room for a significant 'dark' accretion 1 ? In par- 
ticular, can BH merging be the main process building the 
present day SMBH MF? 

• how do the X-ray and optical views of BH mass build- 
up compare? 

• what is the typical Eddington ratio of AGNs when they 
accrete most of their mass? 

• how long does the visible phase of the AGNs last? 

• do we really need a mass-to-radiation conversion effi- 
ciency higher than the standard value of e ~ 0.1, as recently 
argued (see e.g. Elvis, Risaliti & Zamorani 2002)? 

The local SMBH MF, including the contribution from 
the spheroidal components of late-type galaxies, is re- 
estimated exploiting and extending the technique outlined 
by Salucci et al. (1999) and recently presented by Shankar 
et al. (2003). The accreted MF will be derived using up-to- 
date LFs as a function of cosmic time for optically and X-ray 
selected AGNs. 

The paper is organized as follows. In Section 2 we criti- 
cally discuss the relationships among luminosity, mass, and 
velocity dispersion of the spheroidal components of galax- 
ies (L sp h, M sph , and a), and Mbh. In Section 3 we present 
and discuss two estimates (derived via the velocity disper- 
sion distribution function (VDF) and the LF, respectively) 
of the local SMBH MF. In Section 4 the accreted mass den- 
sity is estimated. In Section 5 the accreted mass function 
(AMF) is computed and compared with the local MF of 



1 By 'dark' accretion we mean accretion not traced by either 
optical or X-ray surveys. This is quite different from the often 
used term of "obscured" accretion, which is referred to accretion 
on type 2 AGNs (see e.g. Fabian 2003). 



SMBHs. The accretion history and the AGN visibility times 
are analyzed in Section 6. A discussion of the results and 
the conclusions are presented in Section 7. 

The standard flat ACDM cosmological model has been 
used, with Ho = 70kms _1 Mpc -1 , fi m = 0.3, tt A = 0.7. 



2 CORRELATIONS AMONG SMBH MASS, 
GALAXY LUMINOSITY AND VELOCITY 
DISPERSION 

The SMBHs MF can be derived coupling the statistical in- 
formation on local LFs of galaxies with relationships among 
luminosity (or related quantities, such as stellar mass and 
velocity dispersion) and the central BH mass (see e.g. Salucci 
et al. 1999). 

Since the BH mass correlates with luminosity and ve- 
locity dispersion of the bulge stellar population, we need 
separate LFs for different morphological types (which have 
different bulge to total luminosity ratios), and it is conve- 
nient to use galaxy LFs derived in bands as red as possible, 
where the old bulge stellar populations are more prominent. 

2.1 Bulge luminosity versus black hole mass 

McLure & Dunlop (2002) analyzing a sample of 72 active 
and 20 inactive galaxies found that the central BH mass 
and the total R-band magnitude, Mr, of the bulge are 
strictly related. In particular, considering only inactive el- 
liptical galaxies with accurate measurements of BH mass, 
the relation, converted to Ho = 70, reads 

log(--p-) = -0.50(±0.05)M fl - 2.69(±1.04) , (1) 

Mq 

with a scatter of Alog(MaH) = 0.33. It is worth noticing 
that the relationship has been derived using B-band magni- 
tudes, translated to R-band assuming an average color (B- 
R)=1.57. A larger scatter Alog(MeH) — 0.45 was found by 
Kormendy & Gebhardt (2001), who used a sample including 
also lenticular and spiral galaxies. 

For galaxies observed with a spatial resolution high 
enough to resolve the BH sphere of influence, Marconi & 
Hunt (2003) report a tight relation between the SMBH mass 
and the host galaxy bulge K-band luminosity 

l°g(^F-) = 1.13(±0.12)[log(-^)-10.9]+8.21(±0.07)(2) 

Mq L Kq 

with a scatter A log Mbh = 0.31. Translating Eq. (J to 
the K-band using the colour R — K = 2.6 (Kochanek et al. 
2001, with K — K20 = —0.2) as discussed in Sect. 3.1, it is 
apparent that the Marconi & Hunt (2003) relationship yields 
higher BH masses at fixed luminosity. Correspondingly, the 
derived SMBH mass density is up to a factor of 2 higher 
than obtained from Eq. Q. A closer analysis shows that 
most of the discrepancy is due to SMBH in spiral galaxies 
and can be ascribed to the uncertainty in the evaluation of 
their bulge component. Since most of the local mass density 
is contributed by BHs in E and SO galaxies, we decided to 
exploit the relationship reported in Eq. Q. 

As pointed out by McLure & Dunlop (2002) , their Mr - 
log Mbh relation is compatible with a linear relation between 
BH and spheroidal mass, M sp h. Indeed, inserting the result 
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found by Borriello et al. (2003), M sph /L fl oc L° R 21±0 03 , we 
getMBHOcM^r 12 - 



2.2 Velocity dispersion versus BH mass 

While the presence of a strong correlation between BH mass 
and velocity dispersion of the stellar spheroid, Mbh — o, 
is undisputed (Ferrarese & Merritt 2000; Gebhardt et al. 
2000), the value of its slope is still debated. A detailed anal- 
ysis of the available data by Tremaine et al. (2002) yields 

log(M BH |^-) = 4.02(±0.32) log(o-2oo) + 8.13(±0.06) , (3) 
Ho 

0200 being the line-of-sight velocity dispersion in units of 
200 km s" 1 . The slope is in reasonable agreement with the 
findings of Ferrarese (2002), M BH oc cr 4 ' 58±0 ' 52 . It should 
also be noted that the velocity dispersions used by Tremaine 
et al. (2002) refer to a slit aperture 2r e , while those reported 
by Ferrarese refer to r e /8. The scatter around the mean 
relationship is small, A log Mbh = 0.3, possibly consistent 
with pure measurement errors. 

The low mass and low velocity dispersion regime is quite 
difficult to investigate. The analysis of M33 by Gebhardt et 
al. (2001) yields an upper limit on the BH mass (~ I5OOM0) 
more than 10 times below the central value predicted by 
Eq. However the larger upper limit (~ 3OOOM0) claimed 
by Merritt, Ferrarese & Joseph (2001) is consistent with 
the steeper Mbh — a relation found by Ferrarese (2002). 
Filippenko & Ho's (2003) estimate of the mass of the cen- 
tral BH in the least luminous type 1 Seyfert galaxy known, 
NGC 4395, (M BH 10 4 - 10 5 M Q ) is not inconsistent with 
Eq. However it should be also mentioned that the BH 
mass in this case has been estimated using indirect, rather 
than dynamical arguments. The efforts to detect the so 
called intermediate mass BHs (1O 3 M < Mbh £ 10 6 M Q ) 
in galactic centers and therefore to constrain the very low a 
(< 50kms _1 ) end of the Mbh-o" relation have been recently 
reviewed by van der Marel (2003). 

It is worth noticing that the Mbh-c relationship needs 
not to keep a power-law shape down to low Mbh or a 
values. On the contrary, Granato et al. (2004) presented 
a model for the coevolution of QSOs and their spheroidal 
hosts whereby, in the least massive bulges, the BH growth 
is increasingly slowed down by supernova heating of the in- 
terstellar medium as the bulge mass (hence a) decreases. As 
a result, Mbh is expected to fall steeply with decreasing a, 
for log^kms^ 1 )] ^ 2.1. This model also predicts that the 
observed spread around the mean relationship is a natural 
one, deriving mainly from the different virialization redshifts 
of host halos. 



3 THE LOCAL SMBH MASS FUNCTION 

The local SMBH MF can be estimated either from the local 
LF or from the local velocity dispersion function (VDF) of 
spheroidal galaxies and galaxy bulges, through the Mbh~ 
L S ph or the Mbh-o" relation, respectively. Previous studies 
(Yu & Tremaine 2002; Aller & Richstone 2002) have shown 
that the two methods may yield estimates of the local mass 
density of SMBHs differing by a factor of ~ 2. On the other 
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Figure 1. Galaxy luminosity function estimates converted to R- 
band total magnitudes as described in the text. 
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Figure 2. R-band luminosity function estimates for early-type 
galaxies. Data points from Kochanek et al. (2001). 



hand Ferrarese (2002), McLure & Dunlop (2003) and Mar- 
coni et al. (2004) found very good agreement among the 
results of the two methods. 



3.1 Local luminosity functions of spheroids and 
bulges 

The LFs best suited for our purpose are those in red and 
IR bands, which are more directly linked to the mass in 
old stars. Moreover, we need separate LFs for the various 
morphological types with different bulge to total luminos- 
ity ratios. We will use the K-band LF by Kochanek et al. 
(2001), the AVband LF by Cole et al. (2001), the r* band 
LF by Blanton et al. (2001, 2003), Nakamura et al. (2002), 
and Bernardi et al. (2003). To compare LFs defined in differ- 
ent bands we must set up a common definition of the total 
magnitude/luminosity and of average colours. 

Since we are interested in the total luminosity of 
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Figure 3. R-band local luminosity functions estimates for late- 
type galaxies. 



spheroidal components of galaxies, we adopt as total mag- 
nitudes those obtained with a de Vaucouleurs profile. We 
have therefore corrected by —0.2 the surface brightness lim- 
ited magnitudes, K20, of the 2MASS sample and the Pet- 
rosian magnitudes, used by Blanton et al. (2001; 2003) and 
by Nakamura et al. (2002). Both magnitude systems in fact 
are defined for apertures which contain ~ 80% of the to- 
tal flux for the adopted profile. For the Kron magnitudes 
of Cole et al. (2001), in the K„ band, we used a brighten- 
ing of —0.11, which is required to convert them to an r 1//4 
luminosity profile. 

Magnitudes were converted to the R-band using the 
mean colours R — K s — 2.51 and R — r* = —0.11 (Blanton 
et al. 2001). We assume an error of 0.1 mag. on colours and 
include it in our estimate of the final errors on the SMBH 
MF. 

As illustrated by Fig. the different estimates of the 
LF are in very good agreement with each other, except for 
that by Blanton et al. (2003), which is low at bright magni- 
tudes (by a factor ~ 4 at Mr < —24). Indeed the Schechter 
function adopted by the latter authors falls below their own 
data points at M T * -5 log(H /100) = -23 (cfr. their Fig. 5). 
The classification by Kochanek et al. (2001) allows a clear 
cut distinction between early and late type galaxies. A sim- 
ilar classification has also been proposed by Nakamura et 
al. (2002). Figures [5] and [3] show that the agreement is quite 
satisfactory also for early and late types separately, although 
the early-type LF by Bernardi et al. (2003) misses objects 
fainter than Mr ~ —21 because of their velocity dispersion 
criterion (a > 70kms _1 and S/N > 10). 



3.2 From the local luminosity function to the 
SMBH mass function 

Based on Table 1 of Fukugita, Hogan & Peebles (1998), 
to obtain the LF of the spheroidal components we adopt 
the average -R-band bulge-to-total luminosity ratios y = 
0.85 ±0.05 for early-type galaxies and = 0.30 ±0.05 for 
late-type galaxies. 

The SMBH mass function was then computed convolv- 




Log(M bl /M.) 

Figure 4. Local mass function of SMBHs hosted by early-type 
galaxies. The dot-dashed line shows the estimate obtained from 
the r*-band LF (Nakamura et al. 2002) coupled with the Mbh — 
^bulge relation (see text); the dark gray area shows the estimated 
errors. The dashed line and light gray area refer the MF derived 
using the bivariate dispersion velocity distribution (cfr. Fig. 171. 
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Figure 5. Global SMBH mass function (solid line) and its uncer- 
tainties (dashed lines). The three dots-dashed and the dot-dashed 
lines show the contributions from SMBHs hosted by early- and 
late- type galaxies, respectively. 



ing the LF with the MBH-L sp h relation by McLure & Dun- 
lop (2002) [see Eq. JJJ], assuming a Gaussian distribution 
around the mean with a dispersion Alog(MBH) = 0.33 ^q q^ . 
These uncertainties encompass most of the values quoted in 
the recent literature (see Sect. 2.2). The errors on the SMBH 
MF include the overall uncertainties on the LF, on the bulge 
fractions, on the A4bh-£s P 1i relation and its scatter. The un- 
certainties on the LF, including those on colours, contribute 
about 70% of the error budget. 

The SMBH MF in early type galaxies, shown in Fig. 0] 
has been estimated using the LF of Nakamura et al. (2002) , 
converted to R-band, and Eq. The match with the MF 
computed via the VDF and the Mbh-o" relation (see be- 
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Table 1. Local SMBH mass densities 



Method 



p° H (10 5 M Q Mpc- 3 /i2 ) 



Early Type Galaxies 

r*+^BH - ibulgc 
bivariate VDF + (M m -cr) 
Sheth VDF + (M BH -cr) 



3.1 
3.0 



+0.9 
-0.8 
1.0 



2.8 



Late Type Galaxies 



+ (M E 



-'bulge 



1.1 ±0.5 



low) is very good. The corresponding SMBH mass density 
amounts to p% s (E) = 3.1± j x 1O 5 M /Mpc 3 , in excellent 
agreement with the findings ol McLure & Dunlop (2003) and 
Marconi et al. (2004), and 30% higher than the estimate by 
Yu & Tremaine (2002) and Aller & Richstone (2002). 

The MF ol SMBH hosted by spiral bulges was com- 
puted in the same way, using the LF for late-type galax- 
ies by Nakamura et al. (2002). Their local mass density is 
Pbh(Sp) = (1.1 ± 0.5) x 10 5 M Q /Mpc 3 , bringing the over- 
all mass density to p BH = (4.2 ± 1.1) x 10 5 M Q /Mpc 3 . The 
local number density of SMBHs with M B h > 10 7 M Q is 
nsMBH ^ (1.3±0.25) x 10" 2 Mpc -3 . As illustrated by Fig. El 
the main contribution to the global mass density comes from 
the range 2 x 10 7 < Mbh < 1 x 10 9 Af Q , mostly populated 
by SMBH in early-type galaxies, while less massive BHs are 
preferentially hosted in late type objects. 

Our determination is very close to the result by Mar- 
coni et al. (2004), who used a methodology similar to ours. 
As suggested by Aller & Richstone (2002) the MF can be 
well represented by a four parameter function, which for our 
determination (per unit dlogAfaH) takes the form: 



$(M B 



. a + 1 


/M BH \ 


8' 


J exp 







(4) 



with $„ = 7.7(±0.3) ■ 10" 3 Mpc" 3 , M, = 6.4(±1.1) ■ 10 7 
M Q , a = -1.11(±0.02) and /3 = 0.49(±0.02) (Ho = 
70kms _1 Mpc -1 ). The formula holds in the range 10 6 < 
Mbh/Mq sC 5 x 10 9 . 




Figure 6. Velocity dispersion function. The solid line is the esti- 
mate obtained from the Nakamura et al. (2002) LF coupled with 
the bivariate (luminosity, a) distribution derived from the SDSS 
data in the r*-band; its uncertainty region is shown by the dashed 
lines. The three dots-dashed line is the estimate by Sheth et al. 
(2003). 
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3.3 The local velocity dispersion function 

The local VDF can be derived from the local galaxy LF 
exploiting the luminosity-cr relation (Gonzalez et al. 2000; 
Sheth et al. 2003), well established for spheroidal galaxies 
(Faber & Jackson 1976). The analysis of a sample of 86 
nearby E and SO galaxies, yields (de Vaucouleurs & Olson 
1982; Gonzalez et al. 2000): 

M Bt = (-19.71 ±0.08) - (7.7 ±0.7) log CT200 +5 log ft, (5) 

with h = ifo/100kms -1 Mpc -1 . However, data for larger 
samples suggest a steeper relation. Bernardi et al. (2003), us- 
ing about 9000 early type galaxies selected from the Sloan 
Digital Sky Survey (SDSS), found L r * tx a 391 , where a 
refers to a r e /8 aperture. The VDF of the SDSS has been 
actually obtained with a fixed aperture of 1. 5 and then 
converted to the r e /8 aperture following the conversion sug- 
gested by j0rgensen, Franx & Kjaergaard (1995). 

Estimates of the local VDF have been derived by Shi- 
masaku (1993) and Gonzalez et al. (2000) (see Kochanek 



Figure 7. Estimates of the local mass function of SMBH hosted 
by early type galaxies, derived from the velocity dispersion func- 
tions in Fig. [H] coupled with the M^n-cr relation by Tremaine 
et al. (2002). The gray area represents the uncertainties on the 
estimate based on the bivariate (luminosity, a) distribution. 



2001 for a comprehensive review), and more recently, by 
Sheth et al. (2003) who were the first to allow for the dis- 
tribution (assumed Gaussian with a luminosity dependent 
width) of data points around the best fit relationship. 

To make a fuller exploitation of the data by Bernardi et 
al. (2003) we have used them to derive the bivariate distri- 
bution Pij = p(Li, cjj), yielding the fraction of objects in the 
r*-luminosity bin centered at Li and in velocity dispersion 
bin centered at <7j. The 9000 objects in the samples, cov- 
ering an absolute magnitude range — 18 ^ M r * ^ —27 and 
a velocity dispersion range 1.8 ^ log(<r) ^ 2.7 (<r in km/s), 
have been subdivided in bins of width 0.05 both in M r * (19 
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bins) and in log(cr) (170 bins). The VDF is then estimated 
as: 

rc(oj) = } J Pani, (6) 

i 

rii — n(Li) being the r*-band LF for early type galaxies by 
Nakamura et al. (2002; Fig. El . The resulting VDF is shown 
in Fig. [(j]with its errors, computed using the formula for the 
propagation of errors in a multivariate function with inde- 
pendent random errors in each variable. The uncertainties 
are bigger towards the two extremes, where the number of 
sampled objects decreases, and smaller around the knee of 
the function. We have checked that our results are indepen- 
dent of the bin size. 

Using the A"-band LF (Kochanek et al. 2001) converted 
to the r*-band adopting a colour K — r* — —2.73, appropri- 
ate for early-type galaxies (Blanton et al. 2001; Kochanek et 
al. 2001) we find differences in the VDF of at most 0.15 dex. 
From Fig. |S| it is apparent that our estimate is very close 
to that by Sheth et al. (2003), apart for the more rapid de- 
cline at low velocity dispersions, due to the selection criteria 
adopted by Bernardi et al. (2003), as noted above. 

The contribution from late type galaxy bulges to the 
VDF is rather difficult to assess. In fact, the bulge-to-disk 
mass ratios depend more on morphology than on luminosity 
and on rotational velocity. Although for a given morpholog- 
ical type a correlation between the bulge velocity dispersion 
and the maximum rotational velocity may be expected, the 
use of the Tully-Fisher relation (among luminosity and rota- 
tion velocity) to infer the velocity dispersion is rather unsafe 
(see discussions of Sheth et al. 2003 and Ferrarese 2002) . 



3.4 From the VDF to the SMBH MF 

In order to get an estimate of the local mass function of 
SMBHs, the local VDF for early-type galaxies can be con- 
volved with the Mbh-o" relation of Tremaine et al. (2002). 
The SDSS velocity dispersions (Bernardi et al. 2003) given 
for an aperture of r e /8, have been converted to 2r e aperture 
using Eq. (16) of Tremaine et al. (2002), while the Mbh — cr 
relation of Tremaine et al (2002) has been estimated with 
velocity dispersions taken within an aperture corresponding 
to 2r e ■ We assume a Gaussian distribution of BH masses at 
constant a, with a dispersion A = 0.30^q Q3 dex. 

The SMBH MF estimates derived from the VDF ob- 
tained through the bivariate probability distribution and 
from the VDF by Sheth et al. (2003) are shown in Fig. □ 
The shaded area shows the uncertainties on the former esti- 
mate, including the contributions from errors on both the 
VDF and A. Again, the decline for M B h < 10 7 M is 
due to the incompleteness of the SDSS sample at low ve- 
locity dispersions. The integrated mass density of SMBH 
in early-type galaxies is p BH = 2.8 x lO 5 A/0/Mpc 3 or 
p BH = 3.0±o ° 6 x 1O 5 M /Mpc 3 if we use the VDF by Sheth 
et al. (2003) or that obtained through the bivariate proba- 
bility function, respectively. As noted above, the evaluation 
of the contribution of SMBHs hosted by late-type galaxies 
through this method is hampered by the poor knowledge 
of the local VDF for their bulges. Adopting the temptative 
estimate by Sheth et al. (2003) for the late type galaxy con- 
tribution to the VDF, coupled with Eq. with the same 
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Figure 8. Optical AGN LF at high-z. Data from the SDSS (Fan 
et al. 2001) and GOODS survey (Cristiani et al. 2004). The solid 
line shows the Croom et al. (2003) power-law model at z = 4.5. 



scatter A = 0.3, we get p° BH = 1.2 x 10 5 M Q /Mpc 3 , nicely 
consistent with the estimate derived from the LF. 

Wyithe & Loeb (2003) obtained a lower estimate of the 
total mass density mainly because they neglected the scatter 
A of the Mbh — o relationship. 



4 THE ACCRETED MASS DENSITY 

Soltan (1982) showed that the total mass density accumu- 
lated by accretion on BHs powering QSOs can be deduced 
from QSO counts, under quite simple assumptions. If e is 
the mass to radiation conversion efficiency, the bolometric 
luminosity is 

£boi = eAf acc c 2 (7) 
and the mass accretion rate reads 

Mbh = (1 - e)Macc (8) 

The conversion of luminosities measured in a given band to 
bolometric luminosities requires the knowledge of bolomet- 
ric corrections fcboi (see, e.g., Elvis et al 1994), which may 
depend on luminosity and/or redshift. The mass accreted 
up to the present time by all AGNs brighter than L can be 
written as 

Pacc(> L) = ^ 2 J dz— J dL'k ho i(L',z)n(L',z)L'. (9) 

where n(L' , z) is the comoving luminosity function. As noted 
by Soltan (1982), p acc is independent of Ho and of the QSO 
lifetime. 

The most complete AGN surveys are those at X-ray 
(hard and soft), optical, and radio wavelengths. The latter 
selection is however rather inefficient, since only ~ 10% of 
AGNs are radio loud. 

4.1 Mass accreted on optically selected QSOs 

The 2dF survey (Boyle 2000; Croom et al. 2003) has pro- 
vided an accurate determination of the redshift-dependent 
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LF of optically selected AGNs with M B < —22.5 and 
z < 2.2. Croom et al. (2003) showed that the data are con- 
sistent with pure luminosity evolution of the form Lb(z) = 
Ls(0)xl0°- 21z(5 ' 476 - z) (for a ACDM model with Q m =0.3), 
peaking at z v ~ 2.74 and exponentially declining at higher 
redshifts. Although the luminosity function is poorly known 
for z > 2.4, there is strong evidence (see Fan et al. 2001 
and Osmer 2003 for a recent review), for a rapid decrease 
with increasing redshift of the space density of bright QSO 
for 2^3. More recently, very deep X-ray (Barger et al. 
2003) and optical (Cristiani et al. 2004) surveys have pro- 
vided strong constraints on the space density of less lumi- 
nous QSOs at high redshift. As illustrated by Fig. |H| the 
Croom et al. (2003) power-law model provides a sufficiently 
accurate description also of the data at z ^ 4. 

Inserting such model in Eq. ©, and integrating it up to 
2 = 6, we get, for fcjf ol = 11.8, appropriate for Lb = (L v u)b 
with ub = 6.8 x 10 14 Hz (Elvis et al. 1994) and e = 0.1: 

P T C = 1.4 x 1O 5 ^|M /Mpc 3 . (10) 

with objects at z ^ 2.2 contributing = 0.8 x 

10 5 M Q /Mpc 3 . Using the Boyle et al. (2000) LF, which is 
however inconsistent with high redshift data, p°J?' increases 
by 20%. Thus the mass density accreted on BHs powering 
the optical QSO emission is a factor ~ 3 lower than the 
estimated local SMBH mass density. 

The estimate of p^cc is affected by uncertainties on fcjj o] 
and on e. An upper limit of fcjf ol = 16 can be derived from the 
Elvis et al. (1994) sample. On the other hand, recent data 
point to a lower bolometric correction than used in Eq. I1UI . 
For instance, McLure & Dunlop (2003) find fcjf ol ~ 8, and 
Vestergaard (2003) finds fcf ol = 9.7 for higher redshift QSOs. 
On the whole we attribute to fc bo i an uncertainty of about 
30%. It is worth noticing that no dependence of the bolo- 
metric correction on optical luminosity has been reported. 

The efficiency e of conversion of accreted mass into out- 
going photons can be as high as ~ 0.4 for extreme-Kerr BHs. 
On the other hand, no firm lower limit to e can be set; in 
extreme cases a BH can grow without radiating any photon 
at all. However, the low value of p^cc does not necessarily 
imply a low value of e, since an additional important contri- 
bution to the local BH mass density is expected from highly 
absorbed hard X-ray selected AGNs, contributing a large 
fraction of the X-ray background energy density, but only 
marginally represented in optical surveys. 

4.2 Mass accreted on X-ray selected AGNs 

Recent data (see Hasinger 2003 for a review) have shown 
that a large fraction of the hard X-ray background is con- 
tributed by partially or completely covered AGNs in the lu- 
minosity range log(L 2 -iok c v) = 42 — 44 erg/s and at redshifts 
2 = 0.5-1. A comprehensive study of the redshift-dependent 
hard X-ray AGN LF, including the distribution of the ab- 
sorption column density Nh, has been recently carried out 
by Ueda et al. (2003; U03 from now on). In the following we 
will refer to their LDDE model, with the additional fraction 
of AGN with 24 < log(iVH) < 25 required in order to fit 
the XRB with the most recent normalizations (Vecchi et al. 
1999; Barcons et al 2000; Gilli 2003). The additional AGN 
fraction implies an increase of the mass density by 25%. 
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Figure 9. Accreted mass density as a function of redshift. The 
solid and the dot-dashed lines show the increase, with decreasing 
redshift, of the comoving accreted mass density as inferred from 
the epoch-dependent hard X-ray luminosity function by Ueda et 
al. (2003), with a luminosity dependent bolometric correction (see 
text), and from the optical luminosity function by Croom ct al. 
(2003), respectively. 



Unfortunately the available information on the overall 
spectral energy distribution of hard X-ray selected objects 
(and particularly of the faint ones, which are the most rel- 
evant to estimate the low mass end of the MF) is scanty, 
so that estimates of the bolometric corrections are difficult. 
The bolometric correction, k^~^° ~ 32, derived by Elvis et 
al. (1994), refers to optically bright quasars. Evidences for 
an increase of the hard X-ray to optical luminosity ratio, 
Lhx I iopt with decreasing optical luminosity have been re- 
ported by Vignali et al. (2003) and bolometric corrections, 
fc 2 ~j 10 ~ 12 — 18, substantially smaller than the Elvis et al. 
(1994) value, have been estimated at least for a few Seyfert 
galaxies (Fabian 2003). Moreover, in order to match the op- 
tical LF of Boyle et al. (2000) starting from the hard X-ray 
LF, U03 had to assume that L2 kev tx I^sooa m close agree- 
ment with the observational data by Vignali et al. (2003). 

The luminosity dependence of the X-ray to optical lu- 
minosity ratio does not necessarily imply a higher efficiency 
of low luminosity objects in producing X-ray (compared to 
optical) photons. An alternative explanation, borne out by 
evidences of larger covering factors for Seyfert galaxies com- 
pared to QSOs, is stronger dust extinction for lower lumi- 
nosity objects which are less capable of pushing away the 
surrounding medium, as suggested long ago (Cheng et al. 
1983). 

If the optical/UV bolometric correction is independent 
of luminosity, the U03 relationship between UV and X-ray 
luminosity implies: 

Inserting Eq. 1111 in Eq. assuming e/(l — e) = 0.1, 
and integrating over the luminosity and redshift intervals 
(41.5 ^ log(Z/2-iokev) ^ 46.5 and 2^3) investigated by 
Ueda et al. (2003), we find 
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pf£ ~ 4.1 x 10 5 M Q /Mpc 3 . (12) 

If we extrapolate the LF up to z = 6, we get a mass density 
larger by 15%. 

As a consistency check, we have subtracted the contri- 
bution of Type 2 AGNs, following the prescriptions given 
by U03, in order to get the contribution to the local mass 
density of Type 1 objects only. We find: 

p^ c pcl 1.5 x 10 5 M Q /Mpc 3 , (13) 

in close agreement with the result obtained using the op- 
tical LF [Eq. llUll l. The relatively large contribution of 
the optically selected AGNs to the local BH mass density 
(> 30%), despite their small (< 20%) contribution to the 
intensity of the HXRB, reflects their lower X-ray to opti- 
cal luminosity ratio. Since p-2c pcl <x [(1 — e )/ e ]^%o\ ° an( ^ 
Pace [(1 — e )/ e ]&bol; from the agreement between the two 
estimates we can conclude that the uncertainty on k^~^° is 
similar to that on fcjf ol , i.e. ~ 30%. 

The mass accretion history is illustrated by Fig. [^show- 
ing the increase with decreasing redshift of the comoving 
accreted mass density, p acc (z), as inferred from the optical 
(dot-dashed line) and from the hard X-ray (solid line) epoch 
dependent comoving luminosity function n(L,z): 

/Oacc(z) = 1 2 £ / dz '4-j / dL' khoiL' n(L',z), (14) 
ec J z az J Lmin 

where the X-ray (but not the optical) bolometric correction 
is a function of luminosity, as discussed above. As shown by 
Fig. |U] most of the accretion occurs at z > 1.5 for optically 
selected AGNs, and at z < 1.5 for hard X-ray selected AGNs. 

The close correspondence of the accreted mass density 
inferred from the hard X-ray LF with the local SMBH mass 
density p B n ~ 4.2 ± 1.0 x 10 5 M /Mpc 3 (see Table [TJ for 
e/(l — e) = 0.1 shows that there is no much room for really 
"dark" accretion (i.e. for accretion with radiative efficiency 
e -C 0.1), confirming the findings by Salucci et al. (1999) and 
of Marconi et al. (2004), unless the luminous phases of the 
AGNs are characterized by radiative efficiencies much higher 
than the usually adopted value. But even if e is close to the 
maximum allowed values (~ 0.3-0.4; Thorne 1974) the ac- 
creted mass accounts for ^ 25-30% of the local SMBH mass 
density, and one would be left with the problem of account- 
ing for the correlations between Mbh and the bulge mass or 
velocity dispersion which arise naturally as a consequence 
of feedback associated to radiative accretion (Silk & Rees 
1998; Cavaliere et al. 2002; King 2003; Granato et al. 2004). 

A more explicit test of the role of accretion is obviously 
the comparison, presented in the next Section, of the result- 
ing MF with the local SMBH MF. 



5 THE LOCAL ACCRETED MASS FUNCTION 

The AMF can be derived from the AGN LF once a re- 
lationship between luminosity and BH mass is established 
(Chokshi & Turner 1992). Salucci et al. (1999) compared, 
under plausible assumptions, the accreted MF with the lo- 
cal SMBH MF to infer information on the accretion history. 
This point has been recently reexamined by a number of 
authors (e.g. Yu & Tremaine 2002; Aller & Richstone 2002; 
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Log(M bl /M.) 

Figure 10. Local SMBH MF (solid line), including SMBHs 
hosted by both early- and late-type galaxies, with its la uncer- 
tainty (shaded area), compared with the accreted MF (dashed 
line) estimated from the X-ray LF by U03, using a luminos- 
ity dependent bolometric correction. Such estimate is obtained 
by differentiating the integral mass function [Eq. 1181 1 with 
A = L/L Edd = 1. 
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Figure 11. Comparison of the accreted MF (dashed line) com- 
puted as in Fig. II 1)1 but for A given by Eq. 1191 with the local 
SMBH MF (solid line, with l<r uncertainties represented by the 
shaded area). The dot-dashed line shows the accreted MF of op- 
tically selected QSOs (M B < -22.5). 



McLure & Dunlop 2003; Yu & Lu 2004; Marconi et al. 2004), 
who reached different conclusions. 

Let us assume that the local SMBH mass is mostly due 
to radiative accretion and that the accretion rate Mbh is 
proportional to Mbh (see e.g. Small & Blandford 1992; Cav- 
aliere & Vittorini 2002; Marconi et al. 2004), at least during 
the main accretion phases. A recent analysis of SDSS quasars 
suggests that the two quantities are correlated (McLure & 
Dunlop 2003), although with a huge scatter, at least partly 
due to the uncertainties in BH mass estimates. An almost 
constant Mbh/Mbh is also expected, according to the phys- 
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Log(M oh /M e ) 

Figure 12. Cumulative local SMBH MF (dashed line) with its 
1(7 uncertainties, compared with the cumulative MF of optically 
selected QSOs (Mb < —22.5; dot-dashed line) plus X-ray selected 
Type 2 AGNs. 
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Figure 13. Iso-x 2 contours in the e— Ao plane for the match be- 
tween the accreted and the local SMBH MFs. The contours are 
labelled with their value of A% 2 . The projections of the A\ 2 = 1 
and Ax 2 = 3.84 contours on the axis corresponding to a param- 
eter give the 68% and 95% confidence intervals, respectively, for 
such parameter. 



ical model of Granato et al (2004), during the fast growth 
of the SMBHs and up to the bright quasar phase. 

If Mbh/Mbh is constant, the bolometric luminosity 
grows exponentially (as does the BH mass): 



Lboi(t) 



eM acc c 
Ac 2 



t E 



1 - e 
Mbh(U) exp 



Mbhc 
(t - U) 



with e-folding time 



e/ 



(l-e)A' 



(15) 



(16) 



where A is the average ratio L/L^dd, is is the Eddington 
time and ti is the time when the growth starts. The e-folding 
time equals the Salpeter time if A = 1. 

The growth stops when the SMBH reaches its maxi- 
mum mass, set equal to its present-day mass Mg H , i.e. we 
neglect the mass increase during the declining phase of the 
light curve (see Yu & Lu 2004). The maximum bolometric 
luminosity is then: 



ibol,max(A^BH) — A— 



(17) 



Under these assumptions, the local SMBH MF is related to 
the epoch-dependent LF in a given observational band by 
the energy balance equation: 



dt 



dLk^ \Ln(L, t) — 



dM° H n(M° H ) Me - M° H 1 



(18) 



where L = L mBX {M^). The local MF n(Af£ H ) is straight- 
forwardly obtained differentiating Eq. 1181 with respect to 
Mbh- It is worth noticing that the above equation can be 
also derived starting from the continuity equation (see e.g. 
Yu & Lu 2004). 

In Fig. HOI the estimated AMF derived from the epoch- 
dependent X-ray LF by U03, assuming Eddington limited 
accretion (A = 1), is compared with the local SMBH MF (in- 
cluding SMBHs hosted by both early- and late-type galax- 
ies). Although the two curves are rather close to each other, 
their shapes differ. The fact that the assumption of Edding- 
ton limited accretion lead to an AMF exceeding the SMBH 
MF in some mass range shows that it cannot be true for 
all epochs and/or luminosities. Indeed, low-z/low luminosity 
AGNs are known to be radiating well below the Eddington 
limit (Wandel et al. 1999) , and recent estimates suggest that 
quasars are in a sub-Eddington regime up to z ~ 2 (McLure 
& Dunlop 2003; Vestergaard 2003). 

The match between the AMF and the local SMBH MF 
indeed improves significantly (Fig. lllH if we adopt a redshift 
dependent Eddington ratio of the form: 



if z > 3 



\{z) 



A [(l + z)/4] a ifz<3 



(19) 



with Ao = 1 and a = 1.4. The discrepancy at Mbh *5 10 9 
Mq is only marginally significant being at slightly more 
than la level. On the other hand, the generally higher AMF 
estimate derived from the X-ray, compared to that from 
the optical, LF (see Fig. 1111 . reflects the strong luminos- 
ity dependence of the fraction of Type 2 AGNs, which are 
represented in the X-ray, but not in the optical, LF. X- 
ray surveys (see, e.g., Hasinger 2003) have shown that the 
Type 2 fraction increases from ~ 30% at high luminosities 
(L2_io kev ~ 10 44 ergs _1 ) to ~ 70-80% at low luminosities 
(I/2-iokeV ~ 3 x 10 42 ergs -1 ), consistent with the results of 
optical spectroscopic surveys of complete samples of nearby 
galaxies, without pre-selection (Huchra & Burg 1992; Ho 
et al. 1997). As a check, we have computed and plotted in 
Fig. 1121 the cumulative accreted mass density function ob- 
tained summing the contribution of the optically selected 
QSOs to that of Type 2 X-ray selected AGNs; again, the 
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Figure 14. Contributions to the local SMBH MF as a function 
of redshift, for several values of Af£ H . 



agreement with the local SMBH mass density function is 
very good. 

We checked that a dependence of A on luminosity, as 
suggested by Salucci et al. (1999), rather than on redshift, 
yields an equally good fit. 

Requiring that the AMF matches the local SMBH MF 
we obtain constraints on the radiative efficiency and on the 
maximum value of the Eddington ratio [Eq. 1191 1. A min- 
imum x 2 analysis yields e ~ 0.09 (+0.04, —0.03) and Ao — 
0.3 (+0.3, -0.1) (68% confidence errors; see Fig. EJ. The 
constraints on the parameter a ruling the evolution of the 
Eddington ratio [Eq. 1191 1 are rather loose (0.3 ^ a ^ 3.5). 



r vis (M£j 



t e f In 



^max(-A^Bti) 



t e f In 



Mm 



(22) 



where Mbh" is the BH mass when radiative accretion yields 
luminosity L ^ L m i n . In the case of the X-ray LF of U03, 
the minimum 2-10 keV luminosity included at z ~ 0.8, 
where the contribution to the X-ray background peaks, is 
log(L m i n ) ~ 42-42.4. With the bolometric correction given 
by Eq. 1111 . the corresponding minimum BH mass is: 



Ml 



5 x 10 4 



A 10 42 



Mr. 



ergs" 



(23) 



The minimum BH mass contributing to the LF of optically 
selected QSOs for the bolometric correction by Elvis et al. 
(1994) is : 



A,/ mm 
JU BH 



3 x 1 7 . 
A 



-10 



-0.4(22.5+A/" ,ax ) 



Mr. 



(24) 



-22.5 (Boyle et al. 2000; Croom et al. 2003). 
The condition that the AGN visibility timescales deter- 
mined via the LF and via the local SMBH mass function are 
equal: 



dt I dLn(L,t) 



l <:f 



M° 

dM° H n(M° BH )H^), 



(25) 



add further constraints on e and Ao. Using the U03 LF and 
assuming a 30% uncertainty on the bolometric correction 
for X-ray luminosities, we find an allowed range for e pretty 
similar to that following from the match between the AMF 
and the local SMBH mass function (0.06 ^ e ^ 0.13) while 
the constraints on the Eddington ratio are looser (0.5 ^ 
At, s; 2). 



6 ACCRETION HISTORY AND AGN 
VISIBILITY TIMES 

Replacing to with t in Eq. 1181 and differentiating with re- 
spect to z and Mj H we get the contributions to the n(Mg H ) 
from different cosmic epochs, shown for several values of 
M£ H in Fig. 1141 which evidences that mass is accreted ear- 
lier and more rapidly by the more massive BHs. 

The the mean time spent by an AGN at L ^ L, averaged 
over the MF, writes: 



J*°dtf™dLn(L,t) 
/~ dM° H n(M° H ) 



(20) 



As noted by Yu & Tremaine (2002), (r; um ) is independent 
of e. It depends however on Aw and on A since 



- A M BH 2 
L — c . 

Kbol tE 



(21) 



Adopting the U03 LF, we find (n um ) ~ 1.5 x 10 s yr if we 
adopt a redshift-dependent Eddington ratio [Eq. I19H 1 and 
~ 8 x 10 7 yr if we keep A = Ao. 

This is however a lower limit to the time interval r v i s , 
spent at L ^ -Lmin, the minimum luminosity included in the 
LF. If M/Mbh = const, the duration of the visible phase is 
simply: 



7 DISCUSSION AND CONCLUSIONS 

The analysis reported in the first part of the paper shows 
that the local mass function of SMBHs can be rather accu- 
rately assessed. More in detail, the MF of SMBHs hosted in 
early-type galaxies can be obtained exploiting the velocity 
dispersion or luminosity functions of host galaxies, coupled 
with the Mbh-o" or MBH-L sp h relationships, respectively. 
The results obtained in the two ways are in remarkable 
agreement with small uncertainties up to Mbh ^ 10 9 Mq 
(cfr. Fig. 0. The contribution from SMBHs hosted by late- 
type galaxies is more uncertain, and is mostly confined to 
the low mass end of the MF. 

The overall SMBH mass density amounts to p B H — 
(4.2 + 1.1) x 1O 5 M /Mpc 3 , with a contribution from SMBHs 
in late-type galaxies of ~ 25%. This value of p% H is higher 
than those found by Yu & Tremaine (2002), by Aller & 
Richstone (2002) and by McLure & Dunlop (2003) (who 
have not considered the contribution from SMBH resid- 
ing in late type galaxies), but is in excellent agreement 
with the results by Marconi et al. (2004). The local num- 
ber density of the SMBHs more massive than 10 6 M© is 
ti(Mbh > 10 6 Mq) ~ 1.7 x 10~ 2 Mpc~ 3 , which corresponds 
to the number of bulges and spheroids with M sp h > 5 x 10 s 
M e . 

The Soltan (1982) argument applied to the hard X-ray 
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selected AGNs, allowing for a luminosity-dependent bolo- 
metric correction (U03; Fabian 2003), yields, for a mass to 
radiation conversion efficiency e = 0.1, an accreted mass 
density of p££ ~ 4.1 x lO 5 A/ /Mpc 3 , in close agreement 
with the local SMBH mass density, indicating that most of 
the BH masses were accumulated by radiative accretion, as 
previously concluded by Salucci et al. (1999). Optically se- 
lected QSOs account for only ~ 35% of the total SMBH 
mass density. The dominant contribution comes from Type 
2 AGNs, mostly missed by optical surveys. 

Not only the mass density, but also the MFs of the 
SMBHs and of the accreted mass match remarkably well, 
if we allow for a decrease of the Eddington ratio A = 
L/LiEdd with redshift [Eq. 1191 1. as suggested by observa- 
tions (McLure & Dunlop 2003; Vestergaard 2003). Optically 
selected, Type 1 AGNs account for the high mass tail of the 
AMF, while Type 2 AGNs take over at lower masses, re- 
flecting the strong increase with decreasing luminosity of the 
Type 2 to Type 1 ratio demonstrated by hard X-ray surveys 
(see, e.g., Hasinger 2003) and consistent with the outcome of 
spectroscopic surveys of complete samples of nearby galax- 
ies, without pre-selection (Huchra & Burg 1992; Ho et al. 
1997). 

The alternative possibility that most of the mass has 
been accumulated by "dark" accretion (i.e. accretion unde- 
tectable by either optical or hard X-ray surveys, as in the 
case of BH coalescence), is severely constrained by the above 
results. In order to make room for this possibility one has 
to assume that the radiative efficiency during the visible 
AGN phases is at the theoretical maximum of e ~ 0.3-0.4. 
But even in this case (unless the bolometric correction is 
far lower than currently estimated) the contribution of ra- 
diative accretion to the local SMBH mass density is Z 25%, 
and one is left with the problem of fine tuning the radia- 
tive and non-radiative contributions in order not to break 
down the match with the local SMBH MF obtained with ra- 
diative accretion alone. One would also face the problem of 
accounting for the tight relationships between BH mass and 
mass or velocity dispersion of the spheroidal host, naturally 
explained by feedback associated to radiative accretion. For 
these reasons it is unlikely that the present day SMBH mass 
function has been built mostly through 'dark' accretion or 
merging of BHs. 

If indeed the SMBH MF has to be accounted for by 
radiative accretion, the requirement that it fits together 
with the AMF constrains the radiative efficiency and the 
maximum Eddington ratio to e ~ 0.09 (+0.04, —0.03) and 
A ~ 0.3 (+0.3, -0.1) (68% confidence errors). The condi- 
tion that the mean AGN visibility timescale computed via 
the LF and via the local MF are equal [Eq. 1251 1 yields an 
allowed range for e very close to the above (0.06 ^ e ^ 0.13) 
and looser (but consistent) constraints on Ao- 

The analysis of the accretion history highlights that the 
most massive BHs (associated to bright optical QSOs) ac- 
creted their mass faster and at higher redshifts (typically at 
2 > 1.5), while the lower mass BHs responsible for most of 
the hard X-ray background have mostly grown at z < 1.5 
(see Figs. 151 and !141 . The different evolutionary behaviour of 
the two AGN populations can be understood if gas accretion 
is regulated by star formation and by feedback both from 
supernova explosions and nuclear activity (e.g. Kawakatu & 
Umemura 2003; Granato et al. 2004). In this framework it 



is expected that the hosts of the most massive BHs have the 
oldest stellar populations (Cattaneo & Bernardi 2003). 

The mass-weighted duration of the luminous AGN 
phase is found to be (ri um ) ~ 0.5-1.5X 10 8 yr. Yu & Tremaine 
(2002), using a similar method, got slightly lower values, be- 
cause they used in Eq. 12011 the optical QSO LF and the 
local MF of SMBHs in early-type galaxies, which, as shown 
in Sect. 3, is only partly accounted for by optically selected 
QSOs. Yu & Lu (2004), from a detailed modelling the lumi- 
nosity evolution of individual QSOs, derived a lower limit 
"n um > 4 x 10 7 yr. 

If the accretion rate per unit BH mass, Mbh/MbHj is 
constant (as in the case of Eddington limited accretion) 
during the main accretion phases, the visibility times in- 
crease with BH mass [Eq. 1221 1. consistent with the finding 
by McLure & Dunlop (2003) that the ratio of SMBH in 
their optically selected sample at z ~ 2 to the correspond- 
ing number density at the present day increases with BH 
mass. These authors estimate, for the most massive BHs 
(Mbh 5^ 10 9 ' 5 M©), a lower limit to r v i s (lifetime in their 
terminology) of 10 8 yr. 

Setting A = 1 and e/(l — e) = 0.1 and inserting in 
Eq. the value of Mbh given by Eq. with Mb = 
-22.5, we get r vis (10 9 - 3 M ) ~ 2 x 10 s yr. Inserting instead 
in Eq. (1221 the value of Mbh given by Eq. with log(L) = 
42 we get r vis (10 6 M Q ) ~ 1.2 x 10 8 yr, r v i s (lO 8 M ) ~ 3 x 
10 8 yr and r vis (10 9 - 3 M ) ~ 4.4 x 10 8 yr. Such visibility times 
increase with decreasing redshift for z < 3 [Eq. 1191 1. 

On the contrary, Marconi et al. (2004) infer mean "life- 
times" increasing with decreasing present day BH masses. 
For A = 1 and e = 0.1 they get ~ 1.5 x 10 8 yr for 
Mbh > 1O 9 M and ~ 4.5 x 10 s yr for Af£ H < 10 8 M Q . 
However, the latter "lifetime" corresponds to ~ 11 e-folding 
times, i.e. to a mass increase by a factor 4.5 x 10 4 . But then 
for a large fraction of their "lifetime" low mass BHs are too 
faint to be included in actual LFs [cf. Eq. 1231 and Eq. 1241 1. 

In this context we also note that the amount of time 
spent above the threshold for inclusion in the LF by objects 
reaching large luminosities/masses cannot be constrained by 
the LFs themselves since such objects are too rare to con- 
tribute significantly to the faint end of the LF. 

Our estimates for r v i s are within the broad constraints 
imposed on one side by the QSO evolution timescale (< 
1 Gyr) and on the other side by the proximity effect (> 
10 4 yr; see Martini 2004). The low "lifetimes" inferred (al- 
though with large uncertainties) from clustering properties 
of optical QSOs (Martini & Weinberg 2001; Martini 2004) 
depend basically on the crude assumption that as soon as a 
massive halo virializes, the QSO appears, without any delay. 

Short visibility times (r v i s ^ 0.02 Gyr) are also implied 
by the models of QSO LFs presented by Haehnelt, Natara- 
jan & Rees (1998), and Wyithe & Loeb (2003). However, 
Hosokawa (2002) has shown that, in the same general frame- 
work, the QSO LF can be reproduced with substantially 
higher values of T v j a . Again, short timescales mainly follow 
from the assumption of immediate QSO ignition at the viri- 
alization of the host DM halo. As shown and discussed by 
Monaco et al. (2000) and Granato et al. (2001), a delay be- 
tween virialization and AGN ignition is a key ingredient to 
understand the QSO LF and the relationship between evo- 
lutionary histories of QSOs and of the massive spheroidal 
galaxies hosting them. Granato et al. (2004) have presented 
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a detailed physical model quantifying such delay and its in- 
crease with decreasing mass. 
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